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ABSTRACT 



We report the discovery of photometric oscillations in the host star of the exoplanet WASP-33 b (HD 15082). The data were obtained in 
the R band in both transit and out-of-transit phases from the 0.3-m telescope and the Montcabrer Observatory and the 0.8-m telescope 
at the Montsec Astronomical Observatory. Proper fitting and subsequent removal of the transit signal reveals stellar photometric 
variations with a semi-amplitude of about 1 mmag. The detailed analysis of the periodogram yields a structure of significant signals 
around a frequency of 21 eye d~', which is typical of 6 Scuti-type variable stars. An accurate study of the power spectrum reveals 
a possible commensurability with the planet orbital motion with a factor of 26, but this remains to be confirmed with additional 
time-series data that will permit the identification of the significant frequencies. These findings make WASP-33 the first transiting 
exoplanet host star with 6 Set variability and a very interesting candidate to search for star-planet interactions. 
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1. Introduction 

Over 100 transiting exoplanets have been confirmed to date, 
most of them orbiting solar-type or late-type stars. WASP-33 
b, a gas giant planet showing transits on a fast-rotating main- 
sequence A5 star, represents an uncommon case that offers the 
possibility of studying an intermediate-mass main-sequence star. 
WASP-33 b was first r eported as a transiting planet candidate by 
, IChristian et al.l d2006l) . but it was not officially announced as an 
^ ■ exoplanet until the study of ICollier Cameron et al.l (|2010|). The 
CT) ' relatively long time lapse may be explained because the host star 
!>■ ! (HD 15082, V = 8.3) is a fast rotator (v sint = 86kms~') and 
^ — I ■ this ha mpers precise radial velocity work. Colli er Cameron et al.l 
^ ! (l2O10h carried out a detailed study considering both photometry 
t"^ ■ and spectral line profile variations during transits and established 
t— I ' an upper mass limit of 4.1 Mj for the planet. In addition, the au- 
CD thors presented evidence of non-radial pulsations in the star and 
1 ■ suggested y Dor-type variability. Furthermore, they also found 
| that WASP-33 b orbits the star in retrograde motion and that the 
• i-H , orbit is inclined relative to the stellar equator. 

Here we use observations taken at the amateur Montcabrer 
H \ Observatory and the professional fully robotic Mo ntsec 
. . , Astronomical Observatory - OAdM (IColome et all 120081) . as 
well as additional observations from the Exoplanet Transit 
Database (hereafter ETD, http://var2.astro.cz/ETD). These al- 
low us to provide the first evidence for photometric oscillations 
on the star WASP-33, and to analyse their amplitude and peri- 
odicity. The presence of a large planet close to a star may cause 
tidal effects that are responsible for multiperiodic non-radial pul- 
sations, and in special cases radial pulsations, on its host star 
(ISchuhll2010h . However, there are only a few known exoplanets 
orbiting pulsating stars, such as V391 Pegasi (sdB ty pe), whose 
plane t was discovered u sing the timing method (ISilvotti et al.l 
l2007l) . lBazot et al.1(l2005l) performed the first asteroseismic anal- 
ysis for a solar-like planet host star, A rae, modelling 43 p - 
modes which were previously identified bv lBouchv et al. (2005). 



Recently, Kepler data were used to characterize the exoplanet 
host HAT-P-7 through an ana lysis of its simultaneously discov - 
ered solar-like oscillations (Chris tensen-Dalsgaard et al1 l2010). 
WASP-33 is the first case where 5 Set pulsations have been ob- 
served in a known transiting exoplanet host star. 



2. Observations and photometry 

The first observations in our dataset were obtained from 
Montcabrer Observatory on 2010 August 26, September 7 and 
14 and October 20, using a 0.3-m Schmidt-Cassegrain telescope 
and a SBIG ST8-XME camera with an AO-8 adaptive optics sys- 
tem, working at a 1.03 "per pixel scale. This is an amateur pri- 
vate observatory with the Minor Planet Center code 213, which 
is located in the suburbs of Barcelona, Spain. Photometry on 
two additional nights with very good photometric conditions 
(September 21 and 28, 2010) was obtained from OAdM and 
with a fully robotic 80-cm Ritchey-Chretien telescope and a FLI 
PL4240 2kx2k camera with a plate scale of 0.36"per pixel. All 
photometry described above was carried out in the Johnson R 
band and by defocusing the images to increase the photomet- 
ric precision. Aperture photometry was performed on the im- 
ages, using three or four comparison stars, as available. The 
first transit was observed with the aim of improving the defini- 
tion of the existing light curv es, most of them incomplete so far 
(ICollier C ameron et al. 2010), whereas the follow-up photome- 
try was obtained at several orbital phases in order to further study 
the oscillations that the first photometric dataset seemed to sug- 
gest. 

Additional recent transit photometry of WASP-33 is avail- 
able at the ETD, including some diagrams and a parameter anal- 
ysis of each transit. WASP-33 b transit /?-band photometric data 
from K. Hose (August 23, 2010), F. Hormuth (August 26, 2010) 
and C. Lopresti (September 11, 2010) were used in our analysis. 
All photometry datasets used in this study are listed in Table [T] 
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Table 1. Photometric datasets used in this work. Monte abrer 
Observatory and OAdM data were specially acquired during 
the course of this work. The rest of the photometry is public 
at the ETD. In transit observations, the photometric root mean 
square (rms) residual per measurement and the transit time de- 
viation (TTD) with respect to the mid-transit ephemeris given in 
ICollier Cameron et al.1 d2010) are calculated from the fits repre- 
sented in Fig. [2] 
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Oil 



Author/ 




Transit / 


rms 


TTD 


Observatory 


Date 


Out of tr. 


(mmag) 


(min) 


K. Hose 


23/08/2010 


T 


3.7 


4.10+1.31 


Montcabrer 


26/08/20 10a 


T 


2.8 


9.76 ± 1.01 


F. Hormuth 


26/08/20 10b 


T 


2.7 


3.18 ±0.86 


Montcabrer 


07/09/2010 


OOT 


3.1 




C. Lopresti 


11/09/2010 


T 


3.9 


-1.66+1.62 


Montcabrer 


14/09/2010 


OOT 


2.3 




OAdM 


21/09/2010 


OOT 


1.6 




OAdM 


28/09/2010 


T 


2.1 


2.83 + 0.50 


Montcabrer 


20/10/2010 


T 


1.8 


7.55 + 0.86 



Table 2. Transit parameters fitted using JKTEBOP. The transit 
depth was computed from the fit. The error bars are l-<x uncer- 
tainties. 



Parameter 


Value 


r, + r 2 


0.309 ± 0.007 


r 2 /ri 


0.1066 + 0.0008 


f(°) 


85.8 ± 2.0 


T - 2450000 


5431.8879 + 0.0003 


rms (mmag) 


3.44 


Depth (mmag) 


14.3 ± 0.2 



The root mean square (rms) residuals given in the third column, 
even though they are affected by the pulsations, give an idea of 
the precision of the photometry from each observatory. 



3. Data analysis and discussion 

The transit ph otometry datasets were analysed with the 
JKTEBOP code (ISoufhworfh et alj|2004allbl). which is based o n 
the EBOP code written by Paul B. Etzel dPopper & Etzellfl98ll) . 
Each transit light curve was corrected for a slight trend us- 
ing low-order polynomials. Preliminary fits were run with the 
JKTEBOP by solving for the sum of the radii relative to the or- 
bital semi-major axis (star + planet; r\ +r2), ratio of radii (n/ri), 
inclination, and time of transit. This was done before combining 
all transits by including a shift to normalize their out-of-transit 
level, so that consistent parameters were o btained for all the tran- 
sits. T he period was fixed to that given in lCollier Cameron et alj 
(1201 Ol) . P orb = 1.2198669 d. To properly weigh the observa- 
tions from different authors, we performed a running average 
of the residuals from a preliminary fit in bins of 20 minutes, 
selecting as individual error for each point the standard devi- 
ation of its bin. The best fit to the overall transit photometry 
is shown in Fig. [TJ and the associated parameters are given in 
Table [2] All parameters are consistent wi thin the uncertainties 
with those in lCollier Cameron et al.l J20 10). A pulsation trend is 
clearly present on the residuals of this fit, and as Fig. [2] illus- 
trates, it is already apparent to the eye that all the curves (some 
of them separated by almost 60 d) show "bumps" that always ap- 
pear at the similar orbital phases (note, e.g., the phases just after 
the egress). 
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Fig. 1. Best fit to the overall /?-band transit photometry phase- 
folded (grey symbols) and binned in steps of 0.001 in phase 
(black symbols) for better visualization. The top panel shows the 
residuals of this fit. Error bars are computed as the mean value 
error of each bin. 



To better determine the characteristic period of the observed 
pulsations, we carried out individual fits to each transit by leav- 
ing the inclination and transit time as free parameters and fixing 
the relative radii of the star and the planet according to the global 
result. This was done to minimize the possible effects of transit 
time variations caused by unknown third bodies in the system 
on the residuals and the effect of trend corrections on the transit 
depth. We find that the inclination of each individual fit is well 
within the uncertainty of that determined for the combination of 
the transits. In Table [TJ the transit time deviations of each transit 
are also listed, however, their scarcity and possible systematics 
caused by the pulsation precludes us of drawing any conclusions 
about the presence of additional planets in the system. 

As already shown, the residuals of the transit fits show a clear 
oscillation pattern that is also present in the out-of-transit pho- 
tometry obtained for WASP-33. The presence of oscillations dur- 
ing the entire orbital period, including out-of-transit phases as 
shown in Fig. [3] rules out the possibility of the "bumps" seen in 
transit as being caused by the effect of t he gravity dar kening in a 
rapidly rotating star such as WASP-33 (Barnes 2009). Thus, we 
combined the residuals from the transit fits and the out-of-transit 
light curves in a single curve to better determine the period of the 
oscillations. Figure|4]shows the results of running a periodogram 
on these data with the PERIOD04 code dLenz & Bregerl 12005, 
http://www.univie.ac.at/tops/period04/) that is based on the dis- 
crete Fourier transform algorithm. A main peak around a fre- 
quency of 21 d~' is clearly present in this figure, corresponding 
to a period of about 68 min. However, a close inspection of this 
peak reveals that it is composed of several other peaks that are 
an imprint of the window function associated to our data. The 
second structure around 6 d -1 , which is similar to the typical 
length of the individual photometric series, is probably spurious 
and may be related to the detrending and the flux-shift of the 
different light curves. 

To better analyse the structure of this peak, we fitted the the- 
oretical transit light curve to the original data by adding a si- 
nusoidal modulation of the light of the system, t hus perform- 
ing a phase dispersion minimization algorithm ( Stellingwerf 
1978). The minimization of this parametric fit was performed b y 
means of a Levenberg-Marquardt algorithm dPress et al.H l992). 
Because of the degeneracy on the solutions, we fitted the ampli- 
tude of the modulation and tested a grid of frequency and phase 
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Fig. 2. /?-band transit photometry of WASP-33. The solid line is 
the best fit to each curve with the relative radii of the star and the 
planet fixed to those found in the global fit using the JKTEBOP 
code. The residuals are shown in the upper panels in the same 
order as transits are displayed. See Table[T]for reference. 
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Fig. 3. Out-of-transit photometry of WASP-33 obtained from 
OAdM on 21 September, 2010. 



difference values. A key element to consider is the weighing of 
each data point. At the high precision of our photometry, photon 
noise is not a good measure of the uncertainty of the measure- 
ment because s ystematic effects in the form of correlated noise 
dominate (e.g.. iPont et al.ll2~006l) . Therefore, we applied a two- 
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Fig. 4. Periodogram of the transit residuals and out-of-eclipse 
phases of WASP-33 using PERIOD04 code. Top panel displays 
the window function. 



step process in which we performed the analysis assuming an 
initial constant arbitrary value of the standard deviation <x=l for 
each data point. The best-fitting solution was then used to cal- 
culate a new value for cr as the local value of the rms residual 
(computed as the running average of the residuals in bins of 20 
minutes), and a final fit was run with weights set accordingly. 
As can be seen in Fig. |3J the lowest x 1 is found for a pulsa- 
tion with a semi-amplitude of 0.98±0.05 mmag and a frequency 
of 21.004+0.004 d -1 at a phase difference with the mid-transit 
time of 264+12° (at reference epoch from Table|2]i. This corre- 
sponds to a period of P — 68.56 + 0.02 min. Figure [6] shows the 
overall residuals and the out-of-eclipse photometry phase-folded 
according to this best-fitting period. As a result of the separation 
between observing nights, the main feature in Fig.[5]is composed 
of many equidistant peaks with similar^- 2 . Interestingly, the pe- 
riodogram peak with lower resulting rms differs from the one 
above as it occurs for a pulsation of about 0.86 mmag in semi- 
amplitude with frequency of 21.311+0.004 d corresponding 
to a pulsation period of 67.57 ± 0.02 min. This illustrates the 
relatively strong impact of the weighing scheme because the 
signal is weak and our time series is relatively short. The lat- 
ter frequency is especially interesting as is happens to be com- 
mensurable with the planet's orbital motion with a factor of 

25.997+0.005. 

Line-profile tomography in ICollier Cameron et al.l d2010h 
provided strong evidence for non-radial pulsations in WASP-33 
with a period around one day, similar to those usually found 
in y Dor stars. Moreover, the authors point out the possibil- 
ity that the retrograde orbiting planet could be tidally induc- 
ing them. Both the period of the photo metric oscillations pre- 
sented here and the stellar properties from lCollier Cameron et al.l 
J20T0I) (r eff = 7400 + 200 K, logg = 4. 3 + 0.2) locate WASP- 
33 we ll within the 5 Set instability strip. lHandler & Shobbrookl 
(2002) present a discussion on the different properties of 5 Set 
and y Dor pulsators. Using the formalism there it can be shown 
that the pulsation constant of WASP-33 (log 2wasp-33 = -1.45) 
perfectly corresponds to the 5 Set domain. The power spectrum 
of 6 Set pulsators is usually very ri ch, as illustrated by the 75+ 
frequencies identified for FG Vir (iBreger et al.1 120051) . but as- 
teroseismic modelling is especially difficult for fast rotators as 
WASP-33. Given this evidence, a different scenario consider- 
ing that WASP-33 belongs to the relatively rare class of hybrid 
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Fig. 5. Results from the Levenberg-Marquardt algorithm method 
displaying the structure of the main periodogrampeak. Both am- 
plitude and^- 2 give best fits for 21.00 d _1 , while the rms is mini- 
mized for a pulsation of frequency 21.31 d _1 , which is commen- 
surable with the planet's orbital period. 











1 















-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



Pulsation phase 

Fig. 6. Residuals of the transit fits and out-of-eclipse photometry 
of WASP-33 phase-folded according to the period of the pulsa- 
tion found in this work, i.e., 68.56+0.02 min, (grey symbols) and 
0.01 phase binning (black symbols). The solid curve is the best 
fit sinusoidal modulation with an amplitude of 0.98 mmag. 



pulsators, showing simultaneous 6 Set and y Dor oscillations 
(lHandler et aD2002tlHandbrll2 009). may be possible. 

One scenario that should be explored in spite of all the ev- 
idence is the possibility that the photometric variations are el- 
lipsoidal in nature and are not caused by pulsations. Indeed, the 
tidal bulge travels particularly fast over the stellar surface, be- 
cause the orbital motion of the planet is retrograde with respect 
to the stellar rotation. Note that there is an indeterminacy in the 
stellar rotation velocity because it is not a given assumption that 
the inclination of the stellar spin axis corresponds to the planet's 
orbital inclination. But a simple calculation, equalling the pulsa- 
tion frequency found here to that of the relative orbital motion of 
WASP-33 b above the star surface, renders the ellipsoidal varia- 
tion scenario as physically not valid since the star would have to 
rotate at a largely super-critical speed in terms of gravitational 
break-up. It is more likely that the tidal bulge rotates over the 
stellar surface with a frequency of about 2 d , which is the net 
combination of the orbital and rotation frequencies, and is far 
from the frequencies we find relevant in the periodogram. 



As seen above, the periodogram is still not defined suffi- 
ciently well to assess the true pulsation spectrum of WASP- 
33. However, given the potential interest, we find it appropri- 
ate to briefly address here the possibility of the coupling be- 
tween the pulsation and the orbit. If the pulsation frequency of 
21.311+0.004 d -1 turned out to be real, this could suggest the 
existence of some kind of star-planet interactions that lead to a 
high-order commensurability of 26 between the pulsation and 
orbital periods. Very few cases are known t o date of a S Set 
variable belonging to a close binary system (Wille ms & Claret! 
2005). The best described candidate for tidally-induced oscilla- 
tions is HD 209295, which simultaneously presents y Dor and 
6 Set-type pulsations, and which was photometrically found to 
show several p -modes directly rela ted to the orbital motion of 
its companion ( lHandler et alJBOOl . This poses the tantalizing 
possibility of a similar situation occurring in WASP-33, which 
will certainly be much clearer when additional photometry is 
acquired and puts constraints on the 6 Set frequency spectrum. 

4. Conclusions 

High-precision /?-band photometry has allowed us to present 
the discovery of photometric oscillations in WASP-33, thus be- 
coming the first transiting planet host with 6 Set pulsations. 
These oscillations have a period of about 68 minutes and a semi- 
amplitude of 1 mmag. The frequency spectrum of WASP-33 is 
still not well defined because the size and significance of the 
peaks depends on the scheme used to weigh the data. In a par- 
ticular weighing criterion, the most significant period comes out 
to be commensurable with the orbital period at a factor of 26. 
If this association is assumed to be physically relevant, the mul- 
tiplicity hints at the existence of planet-star interactions. In any 
case, gaining insight into the nature of WASP-33 will necessarily 
require the collection of additional data (possibly multi-colour) 
with longer time baseline and to carry out pulsation and dy- 
namical modelling. The WASP-33 system now represents a new 
benchmark in the world of exoplanets that can provide valuable 
information on stellar pulsations (through, e.g., transit surface 
mapping), on the tidal interactions between planets and stars, 
and on the dynamical evolution of planetary systems. 

Acknowledgements. We thank G. Handler, P. J. Amado and C. Jordi for useful 
discussions. The referee, M. Bazot, is thanked for insightful comments that led to 
the improvement of the paper. Data were partially obtained by J. Martin with the 
Joan Oro Telescope (TJO) of the Montsec Astronomical Observatory (OAdM), 
which is owned by the Consorci del Montsec and operated by the Institute for 
Space Studies of Catalonia (IEEC). This research has made use of the Exoplanet 
Transit Database I http://var.astro.cz/ETD I. 



References 

Barnes, J. W. 2009, ApJ, 705, 683 

Bazot, M., Vauclair, S., Bouchy, E, & Santos, N. C. 2005, A&A, 440, 615 
Bouchy, F. and Bazot, M. and Santos, N. C. and Vauclair, S. and Sosnowska, D. 

2005, A&A, 440, 609 
Breger, M., Lenz, P., Antoci, V., et al. 2005, A&A, 435, 955 
Collier Cameron, A. C, Guenther, E., Smalley, B., et al. 2010, MNRAS, 407, 

507 

Colome, J., Ribas, I., Fernandez, D., et al. 2008, Proc. SPIE 7019, 70192K-1 - 
70192K-12 

Christensen-Dalsgaard, J., Kjeldsen, H., Brown, T. M., et al. 2010, ApJ, 713, 
L164 

Christian, D. J., Pollacco, D. L., Skillen, I., et al. 2006, MNRAS, 372, 1117 
Handler, G. 2009, MNRAS, 398, 1339 

Handler, G, Balona, L. A., Shobbrook, R. R., et al. 2002, MNRAS, 333, 262 
Handler, G. & Shobbrook, R. R. 2002, MNRAS, 333, 251 
Iorio, L. 2010, Ap&SS, 264 I arXiv: 1006.2707 I 

Lenz, P. & Breger, M. 2005, Communications in Asteroseismology, 146, 53 



E. Herrero 1 et al.: WASP-33: The first 6 Scuti exoplanet host 



Pont, R, Zucker, S., & Queloz, D. 2006, MNRAS, 373, 231 
Popper, D. M. & Etzel, P. B. 1981, AJ, 86, 102 

Press, W. H., Teukolsky, S. A., Vetterling, W. T. & Flannery, B. P. 1992, 

Numerical Recipes, Cambridge University Press 
Scargle, J. D. 1982, ApJ, 263, 835 
Schuh, S. 2010, Astronomische Nachrichten, 331, 489 
Silvotti, R., Schuh, S., Janulis, R., et al. 2007, Nature, 449, 189 
Southworth, J., Maxted, P. F. L., & Smalley, B. 2004a, MNRAS, 351, 1277 
Southworth, J., Zucker, S., Maxted, P. F. L., & Smalley, B. 2004b, MNRAS, 355, 

986 

Steffingwerf, R. F. 1978, ApJ, 224, 953 

Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP, 115, 1023 
Willems, B. & Aerts, C. 2002, A&A, 384, 441 

Willems, B. & Claret, A. 2005, Astronomical Society of the Pacific Conference 
Series, 333, 52 



